A historical sequence deletion in a commonly used Bacillus subtilis chromosome integration vector generates undetected loss-of-function mutations.

Since the 1980s, chromosome-integration vectors have been used as a core method of engineering Bacillus subtilis. One of the most frequently used vector backbones contains chromosomally derived regions that direct homologous recombination into the amyE locus. Here, we report a gap in the homology region inherited from the original amyE integration vector, leading to erroneous recombination in a subset of transformants and a loss-of-function mutation in the downstream gene. Internal to the homology arm that spans the 3' portion of amyE and the downstream gene ldh, an unintentional 227 bp deletion generates two crossover events. The major event yields the intended genotype, but the minor event, occurring in ~10 % of colonies, results in a truncation of ldh, which encodes lactate dehydrogenase. Although both types of colonies test positive for amyE disruption by starch plating, the potential defect in fermentative metabolism may be left undetected and confound the results of subsequent experiments.

A historical sequence deletion in a commonly used Bacillus subtilis chromosome integration vector generates undetected loss-of-function mutations.

Since the 1980s, chromosome-integration vectors have been used as a core method of engineering Bacillus subtilis. One of the most frequently used vector backbones contains chromosomally derived regions that direct homologous recombination into the amyE locus. Here, we report a gap in the homology region inherited from the original amyE integration vector, leading to erroneous recombination in a subset of transformants and a loss-of-function mutation in the downstream gene. Internal to the homology arm that spans the 3' portion of amyE and the downstream gene ldh, an unintentional 227-bp deletion generates two crossover events. The major event yields the intended genotype, but the minor event, occurring in ~10% of colonies, results in a truncation of ldh, which encodes lactate dehydrogenase. Although both types of colonies test positive for amyE disruption by starch plating, the potential defect in fermentative metabolism may be left undetected and confound the results of subsequent experiments.

Molecular Time Capsules Enable Transcriptomic Recording in Living Cells.

Live-cell transcriptomic recording can help reveal hidden cellular states that precede phenotypic transformation. Here we demonstrate the use of protein-based encapsulation for preserving samples of cytoplasmic RNAs inside living cells. These molecular time capsules (MTCs) can be induced to create time-stamped transcriptome snapshots, preserve RNAs after cellular transitions, and enable retrospective investigations of gene expression programs that drive distinct developmental trajectories. MTCs also open the possibility to uncover transcriptomes in difficult-to-reach conditions.

BaM-seq and TBaM-seq, highly multiplexed and targeted RNA-seq protocols for rapid, low-cost library generation from bacterial samples.

The ability to profile transcriptomes and characterize global gene expression changes has been greatly enabled by the development of RNA sequencing technologies (RNA-seq). However, the process of generating sequencing-compatible cDNA libraries from RNA samples can be time-consuming and expensive, especially for bacterial mRNAs which lack poly(A)-tails that are often used to streamline this process for eukaryotic samples. Compared to the increasing throughput and decreasing cost of sequencing, library preparation has had limited advances. Here, we describe bacterial-multiplexed-seq (BaM-seq), an approach that enables simple barcoding of many bacterial RNA samples that decreases the time and cost of library preparation. We also present targeted-bacterial-multiplexed-seq (TBaM-seq) that allows for differential expression analysis of specific gene panels with over 100-fold enrichment in read coverage. In addition, we introduce the concept of transcriptome redistribution based on TBaM-seq that dramatically reduces the required sequencing depth while still allowing for quantification of both highly and lowly abundant transcripts. These methods accurately measure gene expression changes with high technical reproducibility and agreement with gold standard, lower throughput approaches. Together, use of these library preparation protocols allows for fast, affordable generation of sequencing libraries.

High-throughput computational discovery of inhibitory protein fragments with AlphaFold.

Peptides can bind to specific sites on larger proteins and thereby function as inhibitors and regulatory elements. Peptide fragments of larger proteins are particularly attractive for achieving these functions due to their inherent potential to form native-like binding interactions. Recently developed experimental approaches allow for high-throughput measurement of protein fragment inhibitory activity in living cells. However, it has thus far not been possible to predict de novo which of the many possible protein fragments bind their protein targets, let alone act as inhibitors. We have developed a computational method, FragFold, that employs AlphaFold to predict protein fragment binding to full-length protein targets in a high-throughput manner. Applying FragFold to thousands of fragments tiling across diverse proteins revealed peaks of predicted binding along each protein sequence. These predictions were compared with experimentally measured peaks of inhibitory activity in E. coli. We establish that our approach is a sensitive predictor of protein fragment function: Evaluating inhibitory fragments derived from known protein-protein interaction interfaces, we found 87% were predicted by FragFold to bind in a native-like mode. Across full protein sequences, 68% of FragFold-predicted binding peaks match experimentally measured inhibitory peaks. This is true even when the underlying inhibitory mechanism is unclear from existing structural data, and we find FragFold is able to predict novel binding modes for inhibitory fragments of unknown structure, explaining previous genetic and biochemical data for these fragments. The success rate of FragFold demonstrates that this computational approach should be broadly applicable for discovering inhibitory protein fragments across proteomes.

Ubiquitous mRNA decay fragments in E. coli redefine the functional transcriptome.

Bacterial mRNAs have short life cycles, in which transcription is rapidly followed by translation and degradation within seconds to minutes. The resulting diversity of mRNA molecules across different life-cycle stages impacts their functionality but has remained unresolved. Here we quantitatively map the 3' status of cellular RNAs in Escherichia coli during steady-state growth and report a large fraction of molecules (median>60%) that are fragments of canonical full-length mRNAs. The majority of RNA fragments are decay intermediates, whereas nascent RNAs contribute to a smaller fraction. Despite the prevalence of decay intermediates in total cellular RNA, these intermediates are underrepresented in the pool of ribosome-associated transcripts and can thus distort quantifications and differential expression analyses for the abundance of full-length, functional mRNAs. The large heterogeneity within mRNA molecules in vivo highlights the importance in discerning functional transcripts and provides a lens for studying the dynamic life cycle of mRNAs.

First-principles model of optimal translation factors stoichiometry.

Enzymatic pathways have evolved uniquely preferred protein expression stoichiometry in living cells, but our ability to predict the optimal abundances from basic properties remains underdeveloped. Here, we report a biophysical, first-principles model of growth optimization for core mRNA translation, a multi-enzyme system that involves proteins with a broadly conserved stoichiometry spanning two orders of magnitude. We show that predictions from maximization of ribosome usage in a parsimonious flux model constrained by proteome allocation agree with the conserved ratios of translation factors. The analytical solutions, without free parameters, provide an interpretable framework for the observed hierarchy of expression levels based on simple biophysical properties, such as diffusion constants and protein sizes. Our results provide an intuitive and quantitative understanding for the construction of a central process of life, as well as a path toward rational design of pathway-specific enzyme expression stoichiometry.

Quantitative Control for Stoichiometric Protein Synthesis.

Bacterial protein synthesis rates have evolved to maintain preferred stoichiometries at striking precision, from the components of protein complexes to constituents of entire pathways. Setting relative protein production rates to be well within a factor of two requires concerted tuning of transcription, RNA turnover, and translation, allowing many potential regulatory strategies to achieve the preferred output. The last decade has seen a greatly expanded capacity for precise interrogation of each step of the central dogma genome-wide. Here, we summarize how these technologies have shaped the current understanding of diverse bacterial regulatory architectures underpinning stoichiometric protein synthesis. We focus on the emerging expanded view of bacterial operons, which encode diverse primary and secondary mRNA structures for tuning protein stoichiometry. Emphasis is placed on how quantitative tuning is achieved. We discuss the challenges and open questions in the application of quantitative, genome-wide methodologies to the problem of precise protein production.

From coarse to fine: the absolute Escherichia coli proteome under diverse growth conditions.

Accurate measurements of cellular protein concentrations are invaluable to quantitative studies of gene expression and physiology in living cells. Here, we developed a versatile mass spectrometric workflow based on data-independent acquisition proteomics (DIA/SWATH) together with a novel protein inference algorithm (xTop). We used this workflow to accurately quantify absolute protein abundances in Escherichia coli for > 2,000 proteins over > 60 growth conditions, including nutrient limitations, non-metabolic stresses, and non-planktonic states. The resulting high-quality dataset of protein mass fractions allowed us to characterize proteome responses from a coarse (groups of related proteins) to a fine (individual) protein level. Hereby, a plethora of novel biological findings could be elucidated, including the generic upregulation of low-abundant proteins under various metabolic limitations, the non-specificity of catabolic enzymes upregulated under carbon limitation, the lack of large-scale proteome reallocation under stress compared to nutrient limitations, as well as surprising strain-dependent effects important for biofilm formation. These results present valuable resources for the systems biology community and can be used for future multi-omics studies of gene regulation and metabolic control in E. coli.

Spurious regulatory connections dictate the expression-fitness landscape of translation factors.

During steady-state cell growth, individual enzymatic fluxes can be directly inferred from growth rate by mass conservation, but the inverse problem remains unsolved. Perturbing the flux and expression of a single enzyme could have pleiotropic effects that may or may not dominate the impact on cell fitness. Here, we quantitatively dissect the molecular and global responses to varied expression of translation termination factors (peptide release factors, RFs) in the bacterium Bacillus subtilis. While endogenous RF expression maximizes proliferation, deviations in expression lead to unexpected distal regulatory responses that dictate fitness reduction. Molecularly, RF depletion causes expression imbalance at specific operons, which activates master regulators and detrimentally overrides the transcriptome. Through these spurious connections, RF abundances are thus entrenched by focal points within the regulatory network, in one case located at a single stop codon. Such regulatory entrenchment suggests that predictive bottom-up models of expression-fitness landscapes will require near-exhaustive characterization of parts.

Sigma factor dependent translational activation in Bacillus subtilis.

Sigma factors are an important class of bacterial transcription factors that lend specificity to RNA polymerases by binding to distinct promoter elements for genes in their regulons. Here we show that activation of the general stress sigma factor, σB, in Bacillus subtilis paradoxically leads to dramatic induction of translation for a subset of its regulon genes. These genes are translationally repressed when transcribed by the housekeeping sigma factor, σA, owing to extended RNA secondary structures as determined in vivo using DMS-MaPseq. Transcription from σB-dependent promoters ablates the secondary structures and activates translation, leading to dual induction. Translation efficiencies between σB- and σA-dependent RNA isoforms can vary by up to 100-fold, which in multiple cases exceeds the magnitude of transcriptional induction. These results highlight the role of long-range RNA folding in modulating translation and demonstrate that a transcription factor can regulate protein synthesis beyond its effects on transcript levels.

Functionally uncoupled transcription-translation in Bacillus subtilis.

Tight coupling of transcription and translation is considered a defining feature of bacterial gene expression1,2. The pioneering ribosome can both physically associate and kinetically coordinate with RNA polymerase (RNAP)3-11, forming a signal-integration hub for co-transcriptional regulation that includes translation-based attenuation12,13 and RNA quality control2. However, it remains unclear whether transcription-translation coupling-together with its broad functional consequences-is indeed a fundamental characteristic of bacteria other than Escherichia coli. Here we show that RNAPs outpace pioneering ribosomes in the Gram-positive model bacterium Bacillus subtilis, and that this 'runaway transcription' creates alternative rules for both global RNA surveillance and translational control of nascent RNA. In particular, uncoupled RNAPs in B. subtilis explain the diminished role of Rho-dependent transcription termination, as well as the prevalence of mRNA leaders that use riboswitches and RNA-binding proteins. More broadly, we identified widespread genomic signatures of runaway transcription in distinct phyla across the bacterial domain. Our results show that coupled RNAP-ribosome movement is not a general hallmark of bacteria. Instead, translation-coupled transcription and runaway transcription constitute two principal modes of gene expression that determine genome-specific regulatory mechanisms in prokaryotes.

Growth-Optimized Aminoacyl-tRNA Synthetase Levels Prevent Maximal tRNA Charging.

Aminoacyl-tRNA synthetases (aaRSs) serve a dual role in charging tRNAs. Their enzymatic activities both provide protein synthesis flux and reduce uncharged tRNA levels. Although uncharged tRNAs can negatively impact bacterial growth, substantial concentrations of tRNAs remain deacylated even under nutrient-rich conditions. Here, we show that tRNA charging in Bacillus subtilis is not maximized due to optimization of aaRS production during rapid growth, which prioritizes demands in protein synthesis over charging levels. The presence of uncharged tRNAs is alleviated by precisely tuned translation kinetics and the stringent response, both insensitive to aaRS overproduction but sharply responsive to underproduction, allowing for just enough aaRS production atop a "fitness cliff." Notably, we find that the stringent response mitigates fitness defects at all aaRS underproduction levels even without external starvation. Thus, adherence to minimal, flux-satisfying protein production drives limited tRNA charging and provides a basis for the sensitivity and setpoints of an integrated growth-control network.

Keeping the Proportions of Protein Complex Components in Check.

How do cells maintain relative proportions of protein complex components? Advances in quantitative, genome-wide measurements have begun to shed light onto the roles of protein synthesis and degradation in establishing the precise proportions in living cells: on the one hand, ribosome profiling studies indicate that proteins are already produced in the correct relative proportions. On the other hand, proteomic studies found that many complexes contain subunits that are made in excess and subsequently degraded. Here, we discuss these seemingly contradictory findings, emerging principles, and remaining open questions. We conclude that establishing precise protein levels involves both coordinated synthesis and post-translational fine-tuning via protein degradation.

Alternative transcription cycle for bacterial RNA polymerase.

RNA polymerases (RNAPs) transcribe genes through a cycle of recruitment to promoter DNA, initiation, elongation, and termination. After termination, RNAP is thought to initiate the next round of transcription by detaching from DNA and rebinding a new promoter. Here we use single-molecule fluorescence microscopy to observe individual RNAP molecules after transcript release at a terminator. Following termination, RNAP almost always remains bound to DNA and sometimes exhibits one-dimensional sliding over thousands of basepairs. Unexpectedly, the DNA-bound RNAP often restarts transcription, usually in reverse direction, thus producing an antisense transcript. Furthermore, we report evidence of this secondary initiation in live cells, using genome-wide RNA sequencing. These findings reveal an alternative transcription cycle that allows RNAP to reinitiate without dissociating from DNA, which is likely to have important implications for gene regulation.

Transcriptional Regulation and Mechanism of SigN (ZpdN), a pBS32-Encoded Sigma Factor in Bacillus subtilis.

Laboratory strains of Bacillus subtilis encode many alternative sigma factors, each dedicated to expressing a unique regulon such as those involved in stress resistance, sporulation, and motility. The ancestral strain of B. subtilis also encodes an additional sigma factor homolog, ZpdN, not found in lab strains due to being encoded on the large, low-copy-number plasmid pBS32, which was lost during domestication. DNA damage triggers pBS32 hyperreplication and cell death in a manner that depends on ZpdN, but how ZpdN mediates these effects is unknown. Here, we show that ZpdN is a bona fide sigma factor that can direct RNA polymerase to transcribe ZpdN-dependent genes, and we rename ZpdN SigN accordingly. Rend-seq (end-enriched transcriptome sequencing) analysis was used to determine the SigN regulon on pBS32, and the 5' ends of transcripts were used to predict the SigN consensus sequence. Finally, we characterize the regulation of SigN itself and show that it is transcribed by at least three promoters: PsigN1 , a strong SigA-dependent LexA-repressed promoter; PsigN2 , a weak SigA-dependent constitutive promoter; and PsigN3 , a SigN-dependent promoter. Thus, in response to DNA damage SigN is derepressed and then experiences positive feedback. How cells die in a pBS32-dependent manner remains unknown, but we predict that death is the product of expressing one or more genes in the SigN regulon.IMPORTANCE Sigma factors are utilized by bacteria to control and regulate gene expression. Some sigma factors are activated during times of stress to ensure the survival of the bacterium. Here, we report the presence of a sigma factor that is encoded on a plasmid that leads to cellular death after DNA damage.

Rapid Accumulation of Motility-Activating Mutations in Resting Liquid Culture of Escherichia coli.

Expression of motility genes is a potentially beneficial but costly process in bacteria. Interestingly, many isolate strains of Escherichia coli possess motility genes but have lost the ability to activate them under conditions in which motility is advantageous, raising the question of how they respond to these situations. Through transcriptome profiling of strains in the E. coli single-gene knockout Keio collection, we noticed drastic upregulation of motility genes in many of the deletion strains compared to levels in their weakly motile parent strain (BW25113). We show that this switch to a motile phenotype is not a direct consequence of the genes deleted but is instead due to a variety of secondary mutations that increase the expression of the major motility regulator, FlhDC. Importantly, we find that this switch can be reproduced by growing poorly motile E. coli strains in nonshaking liquid medium overnight but not in shaking liquid medium. Individual isolates after the nonshaking overnight incubations acquired distinct mutations upstream of the flhDC operon, including different insertion sequence (IS) elements and, to a lesser extent, point mutations. The rapidity with which genetic changes sweep through the populations grown without shaking shows that poorly motile strains can quickly adapt to a motile lifestyle by genetic rewiring.IMPORTANCE The ability to tune gene expression in times of need outside preordained regulatory networks is an essential evolutionary process that allows organisms to survive and compete. Here, we show that upon overnight incubation in liquid medium without shaking, populations of largely nonmotile Escherichia coli bacteria can rapidly accumulate mutants that have constitutive motility. This effect contributes to widespread secondary mutations in the single-gene knockout library, the Keio collection. As a result, 49/71 (69%) of the Keio strains tested exhibited various degrees of motility, whereas their parental strain is poorly motile. These observations highlight the plasticity of gene expression even in the absence of preexisting regulatory programs and should raise awareness of procedures for handling laboratory strains of E. coli.

Production of Protein-Complex Components Is Stoichiometric and Lacks General Feedback Regulation in Eukaryotes.

Constituents of multiprotein complexes are required at well-defined levels relative to each other. However, it remains unknown whether eukaryotic cells typically produce precise amounts of subunits, or instead rely on degradation to mitigate imprecise production. Here, we quantified the production rates of multiprotein complexes in unicellular and multicellular eukaryotes using ribosome profiling. By resolving read-mapping ambiguities, which occur for a large fraction of ribosome footprints and distort quantitation accuracy in eukaryotes, we found that obligate components of multiprotein complexes are produced in proportion to their stoichiometry, indicating that their abundances are already precisely tuned at the synthesis level. By systematically interrogating the impact of gene dosage variations in budding yeast, we found a general lack of negative feedback regulation protecting the normally precise rates of subunit synthesis. These results reveal a core principle of proteome homeostasis and highlight the evolution toward quantitative control at every step in the central dogma.

Genome-Wide Quantitation of Protein Synthesis Rates in Bacteria.

Bacteria produce different amounts of their proteins in response to different conditions. The ability to accurately quantitate the rates of protein synthesis across the genome is an important step toward understanding both underlying regulation and bacterial physiology at a systems level. Ribosome profiling (deep sequencing of ribosome-protected mRNA fragments) enables accurate and high-throughput measurement of such synthesis rates. Ribosomes protect RNAs from nuclease digestion; thus, by collecting and sequencing protected footprints, one can obtain information on the position of every ribosome at the time of cell collection. Assuming ribosomes go on to translate full-length proteins, the density of ribosomes across an ORF can be used to determine protein synthesis rates. Here we outline a step-by-step protocol and discuss the steps where variability and bias may be introduced, including ways to minimize it.